UMBRAL DYNAMICS IN THE NEAR INFRARED CONTINUUM 



A. ANBIC 1 , W. CAO 1 , P. R. GOODE 1 
1. Big Bear Solar Observatory, 40398 North Shore Lane, Big Bear City, CA-92314, USA. 

Submitted to ApJ: 10 December 2010 
Accepted in ApJ: 4 May 2011 

I. Abstract 

We detected peaks of oscillatory power at 3 and ~ 6.5 minutes in the umbra of 
the central sunspot of the active region NOAA AR 10707 in data obtained in the 
near infrared (NIR) continuum at 1565.7 nm. The NIR dataset captured umbral 
dynamics around 50 km below the r 50 o = 1 level. The umbra does not oscillate as a 
whole, but rather in distinct parts that are distributed over the umbral surface. The 
most powerful oscillations, close to a period of ~ 6.5, do not propagate upward. We 
noted a plethora of large umbral dots that persisted for > 30 minutes and stayed in 
the same locations. The peaks of oscillatory power above the detected umbral dots 
are located at 3 and 5 minutes oscillations, but are very weak in comparison with the 
oscillations of ~ 6.5 minutes. 



2. Introduction 

The umbra, the darkest part of a 
sunspot, has complex dynamics. The um- 
bral oscillations, part of the umbral dy- 
namics, have been the subject of many 
research studies. At chromospheric lev- 
els, researchers have observed velocity os- 
cillations with periods close to 3 min- 
utes and velocity amplitudes of about 6 
km/s (Balthasar and Wiehr, 1984; Soltau 
and Wiehr, 1984; Tirospula et al., 2000). 
At the photospheric level, umbral oscil- 
lations have periods closer to 5 minutes, 



with velocity amplitudes of about 75 m/s 
(Lites, 1992). 

From a theoretical standpoint, there 
are two postulated mechanisms for the 
generation of umbral oscillations. One 
is that they are generated on sub- 
photospheric levels by a common ex- 
tended source (Rouppe van der Voort et 
al. 2003; Bogdan and Judge, 2006) and 
propagate upward along the magnetic 
field lines (Lopez Aristre et al., 2001). 
The other postulate suggests that um- 
bral oscillations are generated when slow 
mode waves are trapped within a cavity 
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in the umbra's chromosphere. Broadband 
acoustic noise from the convection zone 
excites the trapped waves (Zhugzhda et 
al.,1983). Lites (1986) found that the 5 
minute photospheric oscillations do not 
drive the 3 minute chromospheric umbral 
oscillations. 

The umbra does not oscillate as the 
whole, but in distinct areas within the 
umbra (Socas-Navarro et al. 2009). 
Socas-Navarro et al. (2009) found that 
the umbra is tremendously dynamic and 
requires a time cadence faster than 20 s to 
resolve the apparent motions of the emis- 
sion source. 

Umbral dots (UDs), small-scale struc- 
tures distributed across the umbra, are 
part of umbral dynamics. Schussler and 
Vogler, (2006) showed that the appear- 
ance of UDs is a natural consequence of 
magnetoconvection under the influence of 
a strong magnetic field. This was con- 
firmed by Cheung et al. (2010) and 
Bharti et al. (2010). Heinemann et al. 

(2007) stated that UDs are caused by 
overturning convection. 

Balthasar and Wiehr (1984) stated 
that the lifetime of UDs is on the or- 
der of 20 minutes. Rimmele (2008), Ri- 
ethmuller et al. (2008), and Ortiz et 
al. (2010) observed the signature of UDs 
dynamics, as described by the model of 
Schussler and Vogler (2006). Rimmele 

(2008) noted that UDs have a lifetime 
close to 30 minutes, as predicted by the 
model of Schussler and Vogler (2006). 
Bharti et al. (2010) similarly predicted 
that the average UD's lifetime is between 
25 and 28 minutes. On the other hand, 



Hamedivafa (2008) stated that the aver- 
age lifetime of the UDs is between 7 and 
10 minutes, while Watanabe et al. (2009) 
measured the average lifetime as 7.3 min- 
utes. 

Socas-Navarro et al. (2009) detected 
a movement of the bright umbral ele- 
ments. They registered horizontal propa- 
gation speeds of 30 km/s and stated that 
their cadence of 20 s is not enough to re- 
solve the fast lateral motion of the oscil- 
latory sources in the umbra. 

In this work, we present observations 
of UDs and umbral oscillations at a level 
of 50 km below the r 50 o = 1- The NIR 
continuum provides easier seeing correc- 
tion with adaptive optics (AO), and much 
lower scattered light, both instrumental 
and atmospheric. Hence we were able 
to observe structures near the telescope's 
diffraction limit for this wavelength even 
without image reconstruction. 



3. Observations and Data analysis 

High-resolution photometric observa- 
tions of the solar active region NOAA 
AR10707 were obtained simultaneously 
in G-band and the near infrared (NIR) 
continuum near 1.6/xm with the Dunn 
Solar Telescope (DST) at the National 
Solar Observatory/Sacramento Peak 
(NSO/SP) on December 1, 2004. Ben- 
efiting from a high order AO system, the 
spatial resolution was close to the diffrac- 
tion limit of the 76 cm aperture DST 
in the NIR continuum, but significantly 
poorer for G-band. 
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To probe the deepest layer of the solar 
atmosphere, exploiting the Sun's opac- 
ity minimum at 1.6//m, NIR photometry 
was employed, which consists of a tune- 
able NIR birefringent Lyot filter devel- 
oped by the Center for Solar- Terrestrial 
Research/NJIT (Cao et al. 2006), an in- 
terference filter with a 5 nm bandpass, 
and a NIR camera. The Lyot filter has a 
very narrow bandpass of 0.25 nm and was 
tuned to the line-free continuum at 1565.7 
nm. The NIR camera (Cao et al. 2005) is 
based on a 1024 x 1024 HgCdTe/A1203 
CMOS focal plane array with a liquid ni- 
trogen cooling system. The pixel size is 
18yumxl8/xm. The output signal is digi- 
talised into 14-bits with a dynamic range 
better than 70 dB. 

The field of view (FOV) was 122" x 122" 
and 105" x 105" for the NIR and the 
G-band, respectively. Ten NIR frames 
were obtained every second, and the best 
frame was selected. The G-band data 
set had a cadence of 0.5 s. All images 
were dark- and flat-field corrected and se- 
lected according to the highest root mean 
square (RMS) contrast. The observa- 
tional run started on 16:54 UT and pro- 
duced datasets 30 minutes long. 

To equalise the time cadence of both 
datasets, we took every second image 
from the G-band dataset, achieving a 1 
s cadence. Thus in both dataset we ac- 
quired a Nyquist period of 2 s. 

The data were co-aligned using a 
Fourier coaligning routine; which uses 
cross-correlation techniques and squared 
mean absolute deviations to provide sub- 
pixel coalignment precision. However, we 



did not implement sub-pixel image shift- 
ing, to avoid substantial interpolation er- 
rors that sometimes accompany the use 
of this technique. The reference image for 
coaligning routine was the floating mean; 
we made a mean image for each interval 
of 30 data images, shifting the interval 
by one image over the time series. De- 
stretching of the images was performed 
to eliminate the influence of seeing dis- 
tortions. The de-stretching routine uses 
bilinear interpolation. The reference im- 
age for the de-stretching routine is the 
floating mean; we made a mean image for 
each interval of 10 data images, shifting 
the interval by one image over the time 
series. 

Wavelet analysis based on work of Tor- 
rence and Compo (1998) was applied. 
The wavelet analysis most useful out- 
put is the power of the detected oscil- 
latory signal that is related to the am- 
plitude of the analysed oscillatory signal. 
(Vigouroux and Delache, 1993; Graps, 
1995; Starch et al. 1997). In our anal- 
ysis we used power of the oscillations as 
the main product of the wavelet analy- 
sis. The oscillatory detection was per- 
formed with use of the Morlet wavelet as 
the mother wavelet: 

(1) Mt)='*-*e iuot e-%, 

where Uq is the non-dimensional fre- 
quency and t the non-dimensional time 
parameter. Morlet wavelet is non- 
orthogonal, based on the Gaussian func- 
tion and very close to the limit of the 
signal processing uncertainty, (y/n). An 
associated Fourier period, P, corresponds 
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to 1.03s, where s is the wavelet scale. 
The wavelet transform is a convolution of 
the time series with the analysing wavelet 
function. The complete wavelet trans- 
form is achieved by varying the wavelet 
scale, which controls both the period and 
temporal extent of the function. Only 
one-dimensional wavelet analysis is done 
for each spatial coordinate. At the begin- 
ning and end of the wavelet transforms 
are regions where spurious power may 
direct result of the finite extent 
of the time series. Those regions are usu- 
ally refereed to as the cone of influence 
(COI), having a temporal extent equal to 
the e-folding time of the wavelet function. 
For the case of our wavelet is: 

p 

(2) t d = V2s = V2 , 

This time scale is the response of the 
wavelet function to noise spikes and is 
used in our detection criteria by requiring 
that oscillations have a duration greater 
than td- The finite nature of our time 
series compelled us to impose the maxi- 
mum period of 12 minutes above which 
we did not accept any detection. The 
approach to automated wavelet analysis 
used here and all the imposed restric- 
tions are detailed in Andic et al. (2010). 
Since the data used contained seeing dis- 
tortions, we took into considerations only 
the oscillations that had 99% confidence 
levels. 

The diameter of the UDs was mea- 
sured as the full-width at half maxi- 
mum (FWHM) of a spatial profile that 
contained the maximum intensity of the 



structure. This analysis was performed at 
the frame with the best contrast. Deter- 
mination of the lifetime and the UDs shift 
in the horizontal plane was performed 
with the nonlinear affine velocity estima- 
tor (NAVE) method (Chae and Sakurai, 
2008), which tracked the detected struc- 
tures through the dataset. 

The areas surrounding the umbra were 
masked out with a binary mask. We 
made an average frame over the time 
series and normalised it to the maxi- 
mum. Then, we multiplied all pixels with 
brightness > 0.75 by and darker pixels 
by 1, masking out in such way everything 
but umbra. 

For the phase analysis of the oscilla- 
tions two methods were applied. One is 
wavelet phase analysis between NIR and 
G-band data sets. This phase analysis 
is similar to the Fourier phase analysis, 
only it provides time and frequency local- 
isation. This localisation is provided by 
the use of wave packets by wavelet anal- 
ysis instead of the infinite wave train of 
Fourier analysis. The difference in cyclic 
phase, A<p can be determined for each fre- 
quency component, u, using the phase in- 
formation contained within the complex 
wavelet transform. The quality of the 
values is represented by phase coherence. 
A time series in this work are extracted 
from the same (x,y) pixel location in both 
datasets. The signals are separated in 
the direction normal to the solar surface. 
Therefore, the phase differences can be 
interpreted as delays caused by the finite 
propagation speed of waves traveling be- 
tween the optical formation heights. The 
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automated method used is described in 
detail in Bloomfield et al. (2004). 

The second method was combining the 
Fourier and Hilbert transform on a single 
dataset (White and Cha, 1973). In short, 
the signal is a complex function. The real 
part is original signal and imaginary part 
is the quadrature of the original signal. 
Since a real function and its quadrature 
are Hilbert transform pairs, the Hilbert 
transform converts one into another. Re- 
sulting transform describes the amplitude 
and phase of a variable in complex plane. 
The signal is transformed into the Fourier 
space then transformed back using the 
Hilbert transform (Stebbins and Goode, 
1987). 

The phase angle spectrum is formed 
in such way that an upward propagat- 
ing wave has a positive phase angle. The 
phase angles are presented in a weighted 
diagram, where weighting is applied per 
sample by cross-power amplitude \fP\Pi 
to produce greyscale phase, Ac/) (Lites 
and Chipman 1979). 

4. Results 

The target of our analysis, the umbra 
of the central spot of AR 10707 (Fig. Q]A.), 
was full of UDs. 

We observed and analysed 70 easily re- 
solved UDs. The UDs were long lived, 
80% of them existed during our entire 
time series of ~ 30 min. During our time 
series, those 80% of UDs stayed in the 
same locations (Figj2]). 



The average size of UDs detected in 
this dataset is 0.86"; double the Dawes' 
limit for the Dunn telescope in NIR con- 
tinuum. Since we analyse objects that are 
close to each other, the Dawes' resolution 
limit is appropriate for our dataset. All 
the sizes of the values below Dawes' limit 
were ignored. 

We analysed the NIR continuum data 
set and G-band data set with wavelet 
analysis. The signal to noise ratio of our 
G-band data was poor since the seeing 
was not good during our observational 
run. Hence we used the G-band data only 
as test data for analysis methods. The 
data sets were 30 minutes long, making 
the highest credible period in this dataset 
12 minutes. With the limitation to the 
period of 12 minutes we avoided the spu- 
rious oscillations that rise in the COI due 
to the finite nature of our time series. To 
avoid noise induced detections near the 
Nyquist frequency, we ignored oscillations 
that did not contain at least 23 points in 
a single period. Thus we analysed oscil- 
lations in the period range of 0.75 to 12 
minutes. 

The NIR continuum shows two strong 
peaks of oscillatory power at 3 and ~ 6.5 
minute (FigfJl solid line). The G-band 
data set (FigJHl dashed line) shows previ- 
ously reported peaks at 3 and 5 minute 
(Balthasar and Wiehr, 1984; Lites, 1992; 
Tirospula et al., 2000) but no significant 
power at ~ 6.5 minute. 
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Figure 1 . Panel A represents the zoomed in part of the umbra, with 
the penumbral elements masked out, from the central spot of the active 
region visible in panel B. Both panels were taken in the NIR continuum. 
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Figure 2. An UD tracked through the time in the NIR dataset. Ar- 
row at each frame marks an example UDs, in this instance twin UDs. 
All selected frames have a contrast above 90% of maximum contrast. 
The marked UDs stay in the same position for the duration of the time 
series. 



The oscillations were separately calcu- 
lated for UDs. The NAVE method pro- 
vides, beside the tracking coordinates of 
the structures, also the intensity changes 
of the structure through the time with the 
intensity error. We performed wavelet 
analysis on each of the UDs intensity 
curves to detect the oscillations emitted 
by UDs (FigHA). The power curve for the 
UDs shows a significantly smaller oscil- 
latory power than the power emitted by 
the whole umbra. The maximum power 
peak for UDs emitted oscillations is at 



~ 1 min, with a broad peak that cov- 
ers periods from ~ 2 to ~ 3 minutes, 
finishing with the peak at ~ 5 min pe- 
riod. The registered power is close to 
the noise level, since the intensity error 
is ~ 20% (FigJU3). The oscillatory curve 
of NIR continuum shows the power peak 
at 3 minutes period with another broad 
peak around 6 minutes. The small peak 
at the 3 minute most probably represents 
oscillations emitted by UD. The other, 
more powerful oscillations registered in 
the other peak might, in part, originate 
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Figure 3. Observed oscillatory power obtained with wavelet analysis, 
both curves normalised at the maximum of the NIR curve. The solid 
line represents the power observed in the NIR continuum, while the 
dashed line represents the power in G-band. 



from a different source. Considering the 
depth of the spectral line formation layer, 
we can speculate that the observed os- 
cillations have a close connection to he- 
lioseismic oscillations. The spectrum of 
the helioseismic oscillations arises from 
modes with periods ranging from about 
1.5 minute to about 20 minute (Gough 
and Toomre, 1991). The observed oscil- 
latory peak may be a consequence of the 
plasma conditions and the shape of the 
cavity below the umbra. 

Due to seeing conditions, the G-band 
data induce too large error into the phase 
difference calculations designed to deter- 
mine if there is wave propagation. The 
umbra observed in our G-band data set 



did not show any distinguishing substruc- 
tures or umbral flashes. We used a com- 
bination of the Fourier and Hilbert trans- 
form to obtain the amplitudes and phases 
of the oscillations on both levels (Fig. EJ. 

At the NIR level, the oscillations with 
significant power have a phase of 0°. The 
range of the phase angle for all registered 
oscillations is from —0.29° to 0°. From 
this range, 89% of pixels in the umbra 
have 0° phase, i.e. the registered waves 
are not propagating in any direction. The 
11% of the pixels in umbra show signals 
of downward propagation with the phase 
angle less than 0.5°, indicating velocities 
< 1 m/s. 
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Figure 4. Panel A presents the oscillatory power calculated for the 
intensity profiles of the individual UDs and averaged over the total 
number of UDs. The power profile shows peaks for 3 and 5 minutes 
oscillations. Panel B presents the averaged intensity variations (blue 
line) of the UDs, with the error bars. The intensity error for UDs 
intensities is 20%. 
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Figure 5. Phase angle of the oscillations registered at NIR contin- 
uum. 

We also calculated the phase differ- detectable wave propagation that might 
ence spectra of the oscillations detected indicate a connection between the oscil- 
in both spectral bands to see is there any lations from both spectral bands. The 
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Figure 6. Phase difference spectrum between oscillations registered 
at G-band level and the NIR level. A small phase angle can be caused 
with the excessive noise in G-band data and not the result of real 
wave propagation. Positive phase difference indicates upward wave 
propagation. 



calculations showed a very slight phase 
difference angle, less than 0.5° (Fig. [H]). 
This small angle might be caused by the 
excessive noise in G-band dataset or by 
the small effects of dissipation through 
the atmosphere and not the true wave 
propagation. Hence we cannot state that 
there is a connection between oscillations 
detected in both spectral bands. 

These findings indicate that there is no 
observable wave propagation between the 
levels. However, due to the high noise lev- 
els in G-band data, this observation has 
to be reevaluated with another dataset. 

The oscillations observed in the NIR 
continuum dataset are distributed over 



the whole umbra. To ensure that the 
oscillations we analysed are real, we 
strengthened the restrictions for auto- 
mated wavelet technique, as well as the 
procedures for the preparation of time 
sets. Use of de-stretching procedure re- 
duces the possibility of intensity oscilla- 
tions induced by shifting the object along 
the X-Y plane (where X and Y stand for 
the spatial coordinates). However, the 
induced intensity oscillations could not 
be 100% removed, because the alignment 
and de-stretching were not performed at 
sub-pixel values. 
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The umbra emits oscillations sporad- 
ically over the time (Fig. [7]). There are 
no constant oscillatory sources. Hence we 
can state that over the time the oscilla- 
tory emission in the NIR continuum level 
is random for almost all periods. How- 
ever, we have to limit this statement. Our 
time series is short and cannot cover the 
longer period oscillations for more than 
few periods; hence we could not draw 
reliable conclusion about the behaviour 
of the longer period oscillations over the 
time. 

We also checked the spatial distribu- 
tion of the dominant oscillation frequency 
across the umbra (Fig. [8]). The um- 
bra does not oscillate whole but 
in patches which are distributed across 
the umbra. For each pixel in the um- 
bra, we marked the period of the regis- 
tered oscillations that carried the most 
power. That period was declared domi- 
nant for that analysed pixel. Each patch 
in the umbra emits oscillations over the 
whole period range we analysed, however 
the maximum power is contained within 
differed period oscillations for the neigh- 
bouring patches. We could not detect a 
connection of the period preference and 
the structures of the umbra. 

The patches are clearly defined and 
separated from each other with a sharp 
change in the oscillatory period. Some 
of the patches that oscillate with the low 
period oscillations have a gradual transi- 
tion from one period oscillations to an- 
other within the patch itself, but still end 
abruptly at the edge. 



Spatially these patches do not corre- 
spond to anything special inside the um- 
bra, and cannot be connected with either 
UD or the spaces between. Some of the 
low period oscillations are connected with 
the observed UDs, however some of the 
UDs have short period oscillations. What 
causes the difference, we do not know. 
We hope that with the higher resolution, 
we will be able to resolve the cause of the 
difference in the patches. 

5. Discussion and Conclusions 

We observed well-defined UDs in the 
NIR continuum level, 50 km below r 50 o = 
1, and we detected oscillations in the um- 
bra with power peaking at 3 and ~ 6.5 
minutes. 

The UDs were distributed across the 
umbra and lasted longer than the dura- 
tion of our dataset (FigfSJ). We did not 
observe the appearance or disappearance 
of existing prominent UDs. Their average 
size is 0.8", and during our dataset they 
tended to stay at approximately the same 
location (Figj2]). Schiissler and Vogler 
(2006) model of magneto-convection in 
the umbra predicts that UDs are caused 
by rising plumes. This aspect was con- 
firmed by the observational work of Rim- 
mele (2008), Riethmuller et al. (2008), 
and Ortiz et al. (2010). In Fig.5 of 
the work (Schiissler and Vogler, 2006), 
we can see the rising of a plume and the 
corresponding appearance of UDs in the 
brightness images. The lifetime of the av- 
erage UD should be around 30 minutes 
of detectable brightness according to this 
model. However from the figure, it is 
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clear that at lower levels in the umbra 
we would be able to see UDs for longer 
time. The example from Schiissler and 
Vogler, (2006) demonstrates that plumes 
exist for almost 40 minutes in the sub- 
photosphere. Moreover, Bharti et al. 
(2010) found that larger UDs live longer. 
We were only capable of resolving large 
UDs, which were stable for the 30 minute 
duration of our data sequence. Thus, we 
can speculate that the longevity of the 
observed UDs is a consequence of efficient 
convection at the NIR continuum forma- 
tion level. Although this line of reasoning 
implies that the same would happen for 
smaller UD's, we could not observe them, 
hence we cannot broaden the statement 
to include the unobserved UDs. 

Our result disagrees with Hamedivafa 
(2008), who found that the lifetime of 
UDs is between 7 and 10 minutes. How- 
ever, Hamedivafa expresses doubt about 
his result, since possible effects were in- 
troduced by the automated analysis used 
for statistical analysis of the UDs. On the 
other hand, Watanabe et al. (2009) also 
used an automated detection code and 
found the mean lifetime to be 7.3 min- 
utes. 

The oscillations observed in the NIR 
continuum have power peaks at 3 and 
~ 6.5 minute period oscillations, while 
the testing of simultaneous observations 
in the G-band show umbral oscillations 
with power peaks near 3 and 5 minutes. 
The observed 3 minute oscillations are 
most probably connected with the os- 
cillations of the UDs (FigHJ), while the 
broader peak around 6 minute oscillations 
might originate from a different source. 



Watanabe et al. (2009) found in UDs 
intensity curves low frequency compo- 
nents close to ~ 8 minutes, while Sobotka 
and Puschmann (2009) reported that 
UDs substructures vary with the time 
scale of ~ 3 minutes. Thus, we can 
speculate that observed power peaks in 
the intensity oscillations could be con- 
nected with same physical mechanism 
that caused the UDs and UD's structures 
to exist and change for the noted pe- 
riod in the spectral lines used in the cited 
works, but at NIR continuum height we 
see those as a pure oscillations and not 
the variations of the UDs and UD's struc- 
tures. 

Oscillations registered in the NIR con- 
tinuum have small negative phase an- 
gle obtained by combining Fourier and 
Hilbert transform. Phase angles close 
to zero indicate evanescent waves. The 
phase angle close to zero is observed for 
the 89% of powerful oscillations indicat- 
ing that oscillations which contribute the 
most to the power peak around 6 minute 
do not propagate upward. The small neg- 
ative phase angle is most probably associ- 
ated with dissipation of evanescent waves. 
We can speculate that these oscillations 
are helioseismic oscillations, which buf- 
fet the photosphere and that these oscil- 
lations do not have a direct connection 
with the upward layers. The typical he- 
lioseismic oscillations have periods rang- 
ing from ~ 1.5 minute to ~ 20 minutes 
(Gough and Toomre, 1991). It is possi- 
ble that the oscillations we observed are 
closely connected with the condition of 
the plasma directly below the umbra. 
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Figure 7. Changes of the observed oscillatory power in the umbra 
over time. We integrated the power for each period individually over 
the whole surface of the umbra. 




Figure 8. Spatial distribution of the dominant oscillation periods 
across the umbra. 
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The phase angle for the most powerful 
oscillations is close to zero and the mini- 
mum oscillatory power is larger than the 
maximum oscillatory power registered in 
G-band. The phase spectra we used in 
this work weights phase angles with the 
power of the oscillatory signal in ques- 
tion, hence phase angles of the weak os- 
cillations will be lost in our figure. We 
observed some waves with the positive 
phase angle, indicating upward propaga- 
tion. However, the weaker power oscilla- 
tions are closer to the sporadic noise de- 
tections of the oscillations and their reli- 
ability is questionable. Thus, we cannot 
state that there is a upward propagation 
of the detected oscillations. 

The peaks at 3 and 5 minutes observed 
in G-band agree with the previous works 
(Balthasar and Wiehr, 1984; Lites, 1992; 
Tirospula et al., 2000), confirming the 
accuracy of our analysis methods. Our 
data set in G-band has poor S/N (sig- 
nal to noise) ratio, thus cannot be used 
for establishing reliable connections be- 
tween NIR dataset findings and the chro- 
mospheric findings. 

The power curve in Fig. [3] for the 
NIR data set shows power for the 3 
minute oscillations small peak. This 
power peak is probably the contribution 
of the oscillations in the UDs. Although 
this might be taken SIS db confirmation of 
the origin of the oscillations from sub- 
photosperic levels (Lopez Aristre et al. 
2001; Rouppe van der Voort et al. 2003; 
Bogdan and Judge, 2006), the low qual- 
ity of the G-band data did not allow us 
to detect the vertical propagation of these 
oscillations and their connection with the 



upper atmospheric levels. Hence we can- 
not present firm proof that umbral oscil- 
lations originate on the NIR continuum 
level. 

The question remains about the other 
possible sources of the registered NIR os- 
cillations. The instrumentation and the 
optical bench did not induce any oscil- 
lations. However, the dataset consist of 
the raw data, and the seeing influences 
differ from visible (G-band) to NIR spec- 
tral lines. Hence, it is necessary to ad- 
just coalignment procedure to include the 
changes in the NIR, as well as in the vis- 
ible part of spectra. Change of the inter- 
val used for the mean reference image for 
the coalignment of dataset caused a shift 
in the oscillatory power peak toward the 
longer periods. With an interval of 10 im- 
ages (i.e. 20 sec), we got a strong power 
peak around the 10 minute period. With 
the interval of 20 images (i.e. 40 sec), 
same peak shifted towards the longer pe- 
riod oscillations peaking at a ~ 13 minute 
period and with the interval of 30 im- 
ages (i.e. 1 min), the peak shifted to 
~ 15 minute period. These oscillations 
were the strongest around the high con- 
trast areas (i.e. an edge of the umbra). 
The position of these oscillations and the 
changes caused by enlarging the image 
interval indicates that seeing distortions 
cause a shift in the images that might 
cause false oscillatory detections. There 
is a possibility that not even with the 30 
images interval used for coalignment, we 
were not completely successful in remov- 
ing all seeing induced oscillations. Thus, 
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these findings have to be reevaluated us- 
ing a dataset where the seeing influences 
are removed. 

Although the whole umbra oscillates, 
the dominant periods of the oscillations 
are not equally distributed over the um- 
bra (Fig. [S]). There are distinct 
over the umbra that oscillate within the 
whole period range we analysed. This re- 
sult agrees, in part, with Socas-Navarro 
et al. (2009), with a slight variation. 
Our patches have different dimensions 
and behave differently than the patches 
described in Socas-Navarro et al. (2009). 
The difference in the oscillatory patches 
in this work and in the work by Socas- 
Navarro et al. (2009) might be caused by 
the different height of formation for the 
spectral lines used. NIR spectral band 
used in this work is located deep in the 
photosphere, and we do not have direct 
connection to the chromospheric obser- 
vations. Thus, we cannot say that these 
patches are same as the ones observed in 
the work by Socas-Navarro et al. (2009). 

The location of the NIR spectral band 
might also lead to the speculation that 
the dominant oscillations with a ~ 6.5 
minute period are helioseismic oscilla- 
tions which do not propagate upward, 
while the oscillations detected in UD that 
have peaks close to 3 and 5 minutes might 
be the oscillations connected with the up- 
per atmospheric layers. However, the 
high error in intensity curves of UD points 
to the necessity of checking this finding 
with some other dataset. 
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